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We study the hydrostatic equilibrium structure of compact stars in the Eddington-inspired Born- 
Infeld gravity recently proposed by Bafiados and Ferreira [Phys. Rev. Lett. 105, 011101 (2010)]. 
We also develop a framework to study the radial perturbations and stability of compact stars in this 
theory. We find that the standard results of stellar stability still hold in this theory. The frequency 
square of the fundamental oscillation mode vanishes for the maximum-mass stellar configuration. 
The dependence of the oscillation mode frequencies on the coupling parameter k of the theory is also 
investigated. We find that the fundamental mode is insensitive to the value of k, while higher-order 
modes depend more strongly on k. 
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I. INTRODUCTION 

General relativity (GR) has been the most successful 
and popular theory of gravity in the past century. From 
its classic predictions on the perihelion advance of Mer- 
cury and deflection of light, to the later predictions such 
as the orbital decay of the Hulse- Taylor binary pulsar due 
to gravitational- wave damping, GR has passed the exper- 
imental observations in these weak-field situations with 
flying colors. Testing the predictions of GR in strong- 
field situations, such as the final stage of binary black 
hole coalescence, will come soon from the detection of 
gravitational waves (see Ref. Q for a review on experi- 
mental tests of GR). 

Despite the great success of GR, the idea that GR may 
not be the correct theory to describe the Universe on cos- 
mological scales is also gaining attention recently. This 
is due to the fact that, if GR is correct, we must then 
require the Universe to be dominated by some unknown 
component, called dark energy, in order to explain the ac- 
celerating expansion of the Universe. In the past decade, 
various alternative theories of gravity which deviate from 
GR on cosmological scales have been proposed in order 
to explain cosmological observations (see Ref. Q for a 
recent review). 

On the other hand, it is also well known that GR is 
not complete because of its prediction of spacetime sin- 
gularities in the big bang and those inside black holes. It 
is generally believed that a consistent theory of quantum 
gravity is needed to resolve this issue. Recently, a new 
Eddington-inspired Born-Infeld (EiBI) theory of gravity 
was proposed by Bahados and Ferreira The appeal- 
ing properties of this theory are that it is equivalent to 
GR in vacuum and can avoid the big bang singularity 
0. The theory differs fr om GR only in the presence of 
matter, and in particular the deviation becomes signif- 
icant at high densities. It is thus reasonable to expect 
that compact stars may be the best astrophysical labo- 
ratories to test EiBI gravity. In fact, Pani et al. [3, [f| 
have recently studied the structure of compact stars in 



EiBI gravity. 

In this work, we report our investigation of compact 
stars in EiBI gravity. We first extend the work of Pani 
et al. 4, 5] by studying the static structure of compact 
stars in this theory using a large set of realistic equations 
of state (EOS). Furthermore, we develop a framework to 
study the radial perturbations and stability of these stars 
in this theory for the first time. 

The plan of the paper is as follows: In Sec. |TT] we 
briefly summarize EiBI gravity. In Sec. IIIII we derive the 
equations for constructing static equilibrium stars in EiBI 
gravity. In Sec. lIVI we present the linearized equations for 
radial oscillations of compact stars. Section [V] presents 
the technique we use to employ realistic EOS models in 
this study. In Sec. IVII we present our numerical results. 
Finally, our conclusions are summarized in Sec. I VIII We 
use units where G = c = 1 unless otherwise noted. 



II. EDDINGTON-INSPIRED BORN-INFELD 
GRAVITY 

Here we briefly summarize the EiBI gravity proposed 
by Bafiados and Ferreira Q . The theory is based on the 
following action 




where l/^l denotes the determinant of / M „ and g = 
R/j,v represents the symmetric part of the Ricci tensor 
and is constructed solely from the connection T^ 7 . The 
matter action Sm is assumed to depend on the metric 
g^v and the matter field ^ m only. Furthermore, it can 
be shown that the action of Eq. ([I} is equivalent to the 
Einstein-Hilbert action when Sm — 0, and hence the 
theory is identical to GR in vacuum 

The constants k and A are related to the cosmological 
constant by A = (A— 1)/k such that, for small nR^ v , the 
action [Eq. (JXJ)] reduces to the Einstein-Hilbert action. 
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We shall fix A = 1 in this work and treat k as the only 
parameter of the theory. The constraints on k based on 
solar observations @, big bang nucleosynthesis and the 
existence of neutron stars 0] have been studied recently. 

In this theory, the spacetime metric g M „ and the con- 
nection are treated as independent fields. The field 
equations are obtained by varying the action [Eq. (JTJ) ] 
with respect to and seperately, and can be writ- 
ten as (with A = 1) Q 



-qqT = 



-gg H 



(2) 
(3) 



where q^ v is an auxiliary metric compatible with the con- 
nection 
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(4) 



and q = |<j^i/|. The stress-energy tensor T^ v satisfies the 
same conservation equations as in GR: 

V„T^ = 0, (5) 

where the covariant derivative refers to the metric g^. 



III. STATIC EQUILIBRIUM 
CONFIGURATIONS 

A. Basic equations 



The structure of compact stars in EiBI theory was first 
studied by Pani et al. [H, [E| . Here we present our formu- 
lation which uses different variables, and hence leads to 
a different set of differential equations which resemble 
more closely the corresponding equations in GR. For a 
static and spherically symmetric spacetime, the space- 
time metric g^ v and the auxiliary metric q^ v are taken 
to be 

g^dx^dx" = -e* (r) di 2 + e x{r) dr 2 + f(r)dn 2 , (6) 
q^daPdsT = -e^dt 2 + e a( - r) dr 2 + r 2 dil 2 . (7) 

Note that the gauge freedom has been used to set qee — 
r 2 . The matter is assumed to be a perfect fluid and is 



described by the stress-energy tensor 

T"- v = {e + P)u^u v + Pg^, 



(8) 



where e and P are the energy density and pressure of 
the fluid, respectively. The four-velocity of the fluid u' 1 
is given by u M = (e~^/ 2 , 0, 0, 0) because of the time and 
spherical symmetry of the spacetime. The stress-energy 
tensor is then simplified to 



T* = -e, T r r = T g e = T£ = P 



(9) 



Note that the indices of T^ v are raised with the spacetime 
metric g^ v . 



From Eq. ([3]), one can obtain the relations 
e p = eH 3 a~ 1 , e a = e x ab. 



(10) 



where a = \fl + Sirne and b = yl — 8ttkP. In addition, 
f(r) is found to be 



r 
ab 



(11) 



The field equations of Eq. ([2]) reduce to the following two 
independent equations 



a 



3 

ab 



2 



(12) 



1 / 1 



k \ab 6 3 ) r 2 r 2 r ^ > ( ) 

where primed quantities denote partial derivatives with 
respect to r. The conservation of the stress-energy tensor 
[Eq. ©] gives 



2P' 

pT 



(14) 



Eqs. (IT2l - (fT4")l can be combined to obtain the following 
two first-order differential equations: 



dP 
dr 



dm 
dr 



1/1 a \ 2m 
2^l^6 + 63^ ) r+ ~ 



1 - 



2m 



P 



1 



1 

4k 



2 - 



3 a 
ab + ¥ 



where the speed of sound c s is calculated by c 2 
The function m(r) is defined by 



/ 2m , 
- _A = 1 )ab. 



(15) 
(16) 
dP/de. 

(17) 



With a given EOS P = P(e), the structure of a hydro- 
static equilbrium star in EiBI gravity is obtained by solv- 
ing Eqs. (|15[) and (|16[) . By expanding a and b in series 
of k, it can be shown that these equations reduce to the 
corresponding structure equations in GR when n — > 
(see, e.g., Ref. 0). 



B. Boundary conditions and numerical scheme 

To construct a static compact star with a given EOS, 
we first integrate Eqs. (fl~5|) and (|16|) by specifying the 
central density e c , and hence the central pressure P c , and 
setting m(0) = 0. The radius of the star R is defined by 
the condition P(R) = 0. As discussed above, EiBI grav- 
ity is equivalent to GR in vacuum. Hence, the interior 
solution should match smoothly to the Schwarzschild so- 
lution at the stellar surface. It can be checked that the re- 
quired conditions are e(R) = P(R) = and hence a(R) = 
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b(R) = 1. We then have e"^ = e- x ^ = (1 - 2M/R), 
where M = m(R) is the mass of the star. We also require 
that = = (1 - 2M/R). The function /9(r), 

and hence ^(r) because of Eq. (|10p . can now be obtained 
by integrating Eq. (fT3|) backward from the surface to the 
center. We have now completed the interior solution of 
the star. It should be noted that the appearance of the 
terms a = VI + 87tk£ and b = y/1 — 8ttkP in the theory 
imposes the following conditions (see also Refs. 0, @): 



We now consider small radial perturbations on the 
static background solution such that F(t,r) — Fo(r) + 
SF(t,r), where F stands for any metric or fluid variable 
and Fq is its background solution. We shall derive the 
linearized field and matter equations by retaining terms 
only of first order in £ and SF. To obtain the oscillation 
mode frequencies, we assume a time dependence e lu>t for 
all the perturbed quantities. Eqs. @ and ([3]) can then 
be reduced to the following equations 



8itkP c < 1 , for k > 0, (18) 
87r|K|e c < 1 , for « < 0. (19) 

IV. RADIAL OSCILLATIONS 
A. Equations for radial perturbations 

In order to check the stability of compact stars con- 
structed in EiBI gravity, we need to study the radial os- 
cillation modes of these stars. The corresponding study 
in GR was first performed by Chandrasekhar almost fifty 
years ago [9| (see also Ref. [l(| for a more recent study) . 
Here we shall derive an eigenvalue equation [see Eq. (|3"0)) 
below] which allows us to obtain the frequencies of radial 
oscillation modes of compact stars in EiBI gravity. Since 
the calculation is somewhat tedious, we shall thus only 
outline the main steps of the derivation in the following. 

We assume that the static background star is per- 
turbed radially so that spherical symmetry is maintained. 
The spacetime metric can still be written as Eq. ^ , but 
now the components gu and g rr depend on both r and t 
(see, e.g., Ref. [HI]): 



g flu dx IJ 'dx h 



f ' r Ut 2 + e x{t ' r) dr 2 + f(r)dn 2 . (20) 



The auxiliary metric g M „ takes the form 

q^dafdx" = -e p{t ' r) dt 2 +e a(t ^dr 2 +2r](t,r)dtdr+r 2 dn 2 , 

(21) 

where the tr component [qt r = rj(t, r)] is in general 
nonzero. The four-velocity of the fluid is now given by 



i" = (-e*/ 2 ,e-* /2 £,0,0), 



(22) 



where £ is the Lagrangian displacement and £ = d£/dt. 
To linear order in £, the nonvanishing components of the 
stress-energy tensor are 



Tf 



TI = Tt 



rpip 
V 



P. 



(23) 
(24) 



5a 

X' 
60 

6p 



-X, 



r 2 Q 2 Se, 



tfi — 4-7TK 

" r Q 2 X 



35P 



Se 



1 - 8irnP 1 + 87TK6 

QaXi 



(25) 
(26) 
(27) 

(28) 



where x = r2 ( e + P)Qi£- in the above, physical quanti- 
ties without "<5" are evaluated on the static background. 
The functions Qi are given explicitly in Appendix [5] On 
the other hand, the linearized conservation equation ([5]) 
becomes 



1 



1 



(e + P)u 2 X = 5P' + -(e + P)5cp l + -(1 



" (29) 

Combining the above linearized equations, we obtain our 
eigenvalue equation for determining the radial oscillation 
modes 



X 



-Wxx ~ W 2 x', 



(30) 



where the functions W\ and W% depend only on the back- 
ground quantities and the frequency square ui 2 (see Ap- 
pendix for their expressions). 



B. Boundary conditions and numerical scheme 

The Lagrangian displacement £ must vanish at the cen- 
ter due to spherical symmetry. From the definition of x, 
this condition is equivalent to 



x(o) = o. 



(31) 



At the stellar surface, the appropriate boundary condi- 
tion is that the Lagrangian variation of the pressure van- 
ishes (SP = 0). It can be shown from Eq. (|26|) . with the 
conditions P(R) — e(R) = 0, that the boundary condi- 
tion is equivalent to the requirement that the displace- 
ment £ is finite at the surface. From the definition of x, 
the boundary condition is thus equivalent to 



where eo and Po refer to the energy density and pressure 
of the unperturbed static background, respectively. It 
should be noted that Eqs. §5§ and imply that the 
auxiliary metric function rj{t, r) is in general nonzero even 
though gtr = 0. 



x(R) = 0. 



(32) 



To find uj 2 numerically, the shooting method is applied. 
We first decompose Eq. (|30p into two first-order differen- 
tial equations for x and x' ■ in practice, we choose a trial 
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eigenvalue u> 2 and start our numerical integration at a 
point near the center. Note that the regularity condition 



14.5 



of x implies that \ 



for small r. We integrate up to 



the stellar surface and check whether Eq. (j52")l is satis- 
fied. The eigenvalue is obtained if the trial ui 2 can satisfy 
the boundary condition. Otherwise, the integration is 
repeated with a different trial to 2 . 



V. EQUATION OF STATE 

In constructing a compact star, one has to specify an 
EOS model that gives the relation between P and e. Al- 
though polytropic EOSs are often used in compact star 
simulations, they are oversimplified and cannot reflect 
the complexity of nuclear matter. On the other hand, 
realistic EOS models are usually presented in tabulated 
forms and one needs in general to perform numerical in- 
terpolations in the study. Alternatively, one can also use 
piecewise polytropic models to fit many tabulated EOSs 
in different density regions [l2j . 

The above two techniques in general work well for 
studying compact stars in GR. However, we found that 
they are not good enough for studying the oscillation 
modes of compact stars in EiBI gravity. The reason 
is that the eigenvalue equation [Eq. ([30)) ] involves the 
derivative c' s , which is proportional to d 2 P/de 2 . It is 
noted that, for the corresponding study in GR, one only 
needs to calculate dP/de ■ In order to apply realistic 
EOS models in our study, we use smooth analytic func- 
tions to model the tabulated EOS models so that the 
second derivative d 2 P/de 2 can be computed analytically, 
and hence numerical errors can be reduced. In particular, 
we use the following analytical representation suggested 
by Haensel and Potekhin fl3l |: 



P = 



a 3 i 3 



f(a 5 (e-a 6 )) 
1 + a^e 

+ (a 7 + a s i)f(a g (a w - e)) 
+ (an + ai 2 e)/(ai 3 (ai4 - e)) 
+ (ai5 + ai 6 e)f(an(ai S - e)), 



(33) 



where P = log(P/dyn cm -2 ), e = log(e/g cm" 3 ), and 
f(x) = l/(e x + l). The 18 constants are fitting param- 
eters. For a given tabulated EOS, we use the Levenberg- 
Marquardt method [l4[ to determine the set of parame- 
ters di that best fits the EOS data points. In this work, 
we consider the following eight realistic EOSs: model 
A [H[, model APR 0, model BBB2 Q3, model C 0, 
model FPS fij^odel SLy4 [H , model UU [Hill, and 
model WS [H|2l|. 

To show the accuracy of our analytical representations 
of the tabulated EOS, we plot in Fig. [Tj the analytical 
fits to the BBB2 (solid line) and FPS (dashed line) EOS 
models. The original EOS data points are denoted by the 
cross (BBB2) and plus (FPS) symbols in the figure. It is 
seen from the figure that the analytic fits match the data 
points of the tabulated EOSs very well. Furthermore, we 




log (e/gcnT) 



FIG. 1: Plot of the analytic fit of the BBB2 and FPS models. 
The crosses and pluses represent the data points in the EOS 
tables and the lines are the analytic fit functions (see also 
Fig. 2 of Ref. [H). 



find that the standard deviations between the original 
EOS data and fitting data points are in general of the 
order 10 -2 or less for all the EOS models. 



VI. RESULTS 

In this section, we shall construct static equilbrium 
compact stars and study their stability in EiBI gravity. 
As discussed in Sec. UU the free parameter of the theory 
is k. For k — 0, the theory is equivalent to GR. On the 
other hand, the more interesting case is k > because 
it is the regime where novel properties of EiBI gravity 
exist. In this regime, the theory leads to a nonsingular 
cosmological model [3j and the existence of pressureless 
stars [3, @ • In this work, we shall consider three different 
values of k defined by 87tk£o = —0.1, 0, 0.1 , where eo = 
I0 15 g cm -3 . These values of k are consistent with the 
recent constraint set by the existence of neutron stars as 
proposed in Ref. 0], though it should be noted that we 
use units where G = c = I in this work. 

To demonstrate the stability of compact stars in EiBI 
gravity, we show in Figs. [2] (a)-2(d) the numerical results 
for the APR, BBB2, FPS, and SLy4 EOS models. For 
each case, we plot the gravitational mass M and the fre- 
quency square lu 2 of the fundamental oscillation mode 
as a function of the central density e c in the upper and 
lower panels, respectively. The circle on each M-e c curve 
corresponds to the maximum- mass stellar configuration. 
It is seen from the figures that the M-e c relation in EiBI 
gravity is qualitatively similar to that in GR. The stellar 
mass increases with the central density until it reaches a 
maximum. In GR, it is known that the mode frequency 
passes through zero at the central density corresponding 
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(a) APR (b) BBB2 




E c ( 10 15 g cm" 3 ) e c ( 10 15 g cm" 3 ) 

FIG. 2: Gravitational mass M and fundamental mode frequency square ui 2 plotted against the central density e c for the four 
EOS models: (a) APR, (b) BBB2, (c) FPS and (d) SLy4. Three different values of K are considered. The circle on each M-e c 
curve corresponds to the maximum-mass configuration. 



to the maximum-mass configuration. This critical den- 
sity corresponds to the onset of dynamical instability in 
the sense that stellar models beyond this critical point are 
unstable against radial perturbations. Our results show 
that this property is also true in EiBI gravity. We see 
that oj 2 also passes through zero at the maximum-mass 
configuration in EiBI gravity. For a given EOS, stellar 
models with central densities less than the critical den- 



sity are stable becuase the fundamental mode frequency 
square w 2 > 0. We have thus demonstrated the stability 
of compact stars in EiBI gravity. 

For a given EOS model, we see that a negative value of 
k in general decreases the maximum mass of a neutron 
star compared to the case of GR. On the other hand, for a 
positive value of k, EiBI gravity can lead to a much larger 
maximum mass. As pointed out in Ref. Q, this has 
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14 -3 

log(e./10 gem ) 



FIG. 3: Upper panel: Gravitational mass as a function of 
the central density for compact stars modeled by the APR 
EOS and 87r«;eo = 0.1. The curve terminates at the point C 
where Eq. (|18|) is violated. Lower panel: Frequency square of 
the fundamental (solid line) and first harmonic (dashed line) 
modes as a function of the central density. 



TABLE I: Mass and radius of the maximum-mass stellar con- 
figuration for each EOS model. Mass M max is expressed in 
solar mass units, and radius R is expressed in km. 





87Tfi;eo = 


= -0.1 


87r«;eo 


= 


87TK60 


= 0.1 


EOS 


Mmax 


R 


Mmax 


R 


Mmax 


R 


A 


1.240 


7.65 


1.656 


8.21 


1.916 


8.98 


APR 


1.922 


9.23 


2.191 


9.82 


2.375 


10.47 


BBB2 


1.585 


8.83 


1.922 


9.36 


2.152 


10.03 


C 


1.488 


9.29 


1.838 


9.68 


2.089 


10.33 


FPS 


1.434 


8.60 


1.800 


9.10 


2.050 


9.79 


SLy4 


1.744 


9.33 


2.052 


9.86 


2.268 


10.51 


uu 


1.932 


9.08 


2.196 


9.67 


2.371 


10.34 


ws 


1.495 


8.90 


1.845 


9.42 


2.094 


10.05 



the interesting implication that some softer EOS models, 
which are ruled out in GR by the recent discovery of a 
neutron star with M w 1.97M© [23j . would be revived 
in EiBI gravity. For example, the maximum mass of a 
neutron star for the FPS EOS in GR is M = 1.8M©. 
However, it can increase to 2M© in EiBI gravity for the 
case 87TK£o =0.1. In Table HI we list the mass and radius 
of the maximum-mass configuration for each EOS model 
we consider in this work. Furthermore, we list the mass 
and fundamental mode frequency as a function of the 
central density for each EOS model in Appendix [B] 

It is recalled that Eqs. (fF8"|) and p^|) must be fulfilled in 
order to construct compact stars in EiBI gravity. In the 
upper panel of Fig. [3] we plot M against e c for compact 
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FIG. 4: Frequency square ui 2 plotted against the number of 
nodes N in the mode eigenfunctions for compact stars mod- 
eled by the APR EOS. The stellar mass M = 1.25M is fixed. 

stars modeled by the APR EOS and 87tk£o = 0.1 again, 
but now the central density is extended to the point C 
where Eq. (|18[) is violated. In the lower panel, we plot 
the frequency square w 2 of the fundamental and first har- 
monic modes against e c . As we have seen in Fig. [2] (a), 
the fundamental mode becomes unstable at the first crit- 
ical point A. However, we now also see that there exists 
a second critical point B at higher densities. Figure [3] 
shows clearly that it is the first harmonic that changes 
stability at this point. The stellar models from B to C 
are still unstable even though dM/de c > in this range. 
Similar to standard neutron stars in GR Q , we have thus 
seen that the criterion dM /de c > also does not guar- 
antee stellar stability in EiBI gravity. 

To end this section, we study how the value of k would 
affect the oscillation mode frequencies. Since the gravi- 
tational mass of a compact star can usually be measured 
more accurately than its radius (see Ref. [24j for a review 
on neutron star observations) , we shall thus consider the 
oscillation modes of compact star models with a given 
mass for different values of k. In Fig. 21 we plot the fre- 
quency square uj 2 against the number of nodes N in the 
mode eigenfunctions. The stellar models have the same 
mass M = 1.25M© and are described by the same APR 
EOS. Three different cases 87tk£o = —0.1, 0, 0.1 are con- 
sidered as before. Note that the fundamental mode has 
zero nodes, and we see from the figure that the frequency 
of this mode is insensitive to the value of k. On the con- 
trary, higher-order modes (i.e., modes with larger values 
of N) depend more sensitively on the value of n. In par- 
ticular, a positive (negative) value of k would decrease 
(increase) the frequencies of higher-order modes. While 
we only show the results of three different values of k in 
Fig. 21 we have seen that this property is in general true 
for other values. 
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Finally, we note that the sensitivity of the frequencies 
of higher-order modes on k may be understood heuristi- 
cally by noting that, in the nonrelativistic limit of EiBI 
gravity, the standard Poisson equation in Newtonian the- 
ory is modified by having an extra source term of the form 
«V 2 p, where p is the mass density A larger variation 
of density perturbation, which is the case for higher-order 
modes, would then lead to a larger difference in the stellar 
dynamics comparing to the Newtonian results. On the 
other hand, the fundamental mode eigenfunction has no 
node and is monotonically increasing in the star. Hence, 
the magnitude of ~S7 2 5p (with 5p being the density per- 
turbation associated to the mode) should in general be 
smaller than those of higher-order modes. This makes 
the frequency of the fundamental mode not so sensitive 
to k. We believe that this property is also true in our 
fully relativistic study, as shown in Fig. [4] 



VII. CONCLUSION 

In this paper, we have studied the equilibrium struc- 
ture of compact stars in EiBI gravity using a set of eight 
realistic EOS models. Our formulation of the structure 
equations is different from the recent works of Pani et 
al. [3, [|| in such a way that the resulting differential equa- 
tions resemble more closely the corresponding equations 
in GR. We solved the structure equations numerically 
and found that the maximum mass of a neutron star can 
be larger than that in GR when the parameter n in EiBI 
gravity is positive. This implies that softer EOS models, 
which are ruled out in GR by the recent discovery of a 
neutron star with mass nearly 2Mq , would be revived in 
EiBI gravity @. 

We have also developed a theory of radial perturbation 
and studied the stability of compact stars in EiBI grav- 
ity by calculating the oscillation mode frequency square 
w 2 for the first time. In contrast to the situation in GR, 
since the oscillation equation in EiBI gravity involves the 
second derivative dP 2 /de 2 , we found that using standard 
techniques such as numerical interpolation or piecewise 
polytropic representation [12[ to handle realistic EOS ta- 
bles in the calculation would not produce reliable numer- 
ical results. This is because the data points in standard 
EOS tables are usually not dense enough to allow an 
accurate calculation of dP 2 /de 2 . We thus followed Ref. 
[HI and used an 18-parameter analytic representation to 
model the EOS in our calculation. It should, however, be 
emphasized that the analytic fitting is not fundamentally 
essential to our stability analysis. It is employed in this 
work in order to improve numerical accuracy only. 

We find that the standard results of stellar stability 
still hold in EiBI gravity. For a sequence of stars mod- 
eled by the same EOS, we found that the fundamental 
mode frequency square passes through zero at the central 
density corresponding to the maximum-mass configura- 
tion. Similar to the analysis of compact stars in GR, this 
point marks the onset of instability. Stellar models with 



central densities less than the critical point are stable 
because uj 2 > 0. Furthermore, we also found that the 
criterion dMj de c > does not guarantee stellar stability. 

We have also studied the effects of the parameter k on 
the oscillation modes. For a fixed stellar mass, we found 
that the fundamental mode frequency is insensitive to the 
value of n. On the contrary, the frequencies of higher- 
order modes depend strongly on k. In particular, a pos- 
itive (negative) value of n would decrease (increase) the 
frequencies of higher-order modes. Our results thus sug- 
gest that the detection of higher-order radial oscillation 
modes might provide useful constraints on the value of k 
if they could be excited to large amplitudes. Of course, 
in reality oscillation modes could in general be excited 
to large amplitudes only in some catastrophic situations 
such as core collapse supernovae. However, those catas- 
trophic events are in general highly nonspherical, and 
hence nonradial oscillation modes could become more rel- 
evant. Studying nonradial oscillations of compact stars 
and the corresponding gravitational wave signals would 
be a natural extension of our present work. We hope 
to return to this issue in the future and study whether 
the gravitational wave signals emitted by compact stars 
could be used to constrain EiBI gravity. 



Appendix A: List of functions 

Here we present the expressions for the various func- 
tions Qi and W% used in Sec. [TV] 
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Appendix B: Numerical data 

In this appendix, we provide the numerical data for 
the gravitational mass M and fundamental mode fre- 
quency square ui 2 . In Tables [TTI to IIX1 we list M and 
uj 2 as a function of the central density e c for three differ- 
ent parameters: 87tk£o = —0.1, 0, 0.1. In the GR limit 
(k = 0), we have checked that stellar models obtained 
by using our analytical representations of the EOS and 
those by the original EOS data (with numerical interpo- 
lation) agree very well. For example, the central densities 
of the maximum-mass configurations obtained by the two 
methods are in general different by about 1% only. 



TABLE III: Model APR. 







- -0.1 


87r«;eo = 


87TKeo = 0.1 




M 


uj 2 


M 


u? 


M 


w 2 


10 15 g cm -3 


Mq 


10 6 s- 2 


M Q 


10 6 s- 2 


Mq 


10 6 s" 2 


2.00 


1.850 


186.4 


2.134 


152.0 


2.351 


97.1 


2.20 


1.882 


139.1 


2.165 


104.6 


2.371 


43.2 


2.40 


1.902 


98.8 


2.182 


63.8 


2.375 


-5.0 


2.50 


1.909 


80.8 


2.187 


45.5 


2.373 


-27.3 


2.60 


1.914 


64.0 


2.190 


28.4 


2.370 


-48.7 


2.70 


1.917 


48.4 


2.191 


12.3 


2.364 


-69.2 


2.80 


1.919 


33.8 


2.192 


-2.7 


2.358 


-88.8 


2.90 


1.921 


20.2 


2.191 


-16.9 


2.350 


-107.8 


3.00 


1.922 


7.4 


2.189 


-30.2 


2.341 


-126.2 


3.10 


1.922 


-4.7 


2.187 


-42.8 


2.331 


-144.1 



TABLE II: Model A. 





87TK60 


= -0.1 


87rfi;eo = 


87rreeo = 0.1 




M 


uj 2 


M 


a; 2 


M 


u, 2 


10 15 g cm" 3 


M e 


10 6 s- 2 


M Q 


10 6 s- 2 


Mq 


10 6 s- 2 


3.30 


1.234 


92.8 


1.639 


112.5 


1.913 


42.9 


3.40 


1.236 


76.3 


1.644 


98.2 


1.915 


24.8 


3.50 


1.238 


60.2 


1.647 


84.3 


1.916 


7.0 


3.60 


1.239 


44.4 


1.650 


70.8 


1.916 


-10.7 


3.70 


1.240 


29.1 


1.652 


57.5 


1.915 


-28.1 


3.80 


1.240 


14.1 


1.654 


44.7 


1.913 


-45.5 


3.90 


1.240 


-0.5 


1.655 


32.1 


1.911 


-62.6 


4.00 


1.240 


-14.7 


1.656 


19.9 


1.909 


-79.6 


4.10 


1.240 


-28.5 


1.656 


7.9 


1.906 


-96.5 


4.20 


1.239 


-42.0 


1.656 


-3.8 


1.902 


-113.3 



TABLE IV: Model BBB2. 





8"KK€o 


= -0.1 


87r/teo = 


87TK60 =0.1 




M 


uj 2 


M 


J 1 


M 


w 2 


10 15 g cm" 3 


Mq 


10 6 s- 2 


Mq 


10 6 s- 2 


Mq 


10 6 s- 2 


2.30 


1.550 


143.4 


1.880 


133.9 


2.130 


93.8 


2.40 


1.559 


125.1 


1.891 


116.6 


2.139 


74.1 


2.50 


1.567 


107.4 


1.900 


99.8 


2.145 


55.0 


2.60 


1.573 


90.3 


1.907 


83.7 


2.149 


36.3 


2.70 


1.577 


73.8 


1.913 


68.1 


2.152 


18.1 


2.80 


1.580 


57.8 


1.917 


53.2 


2.152 


0.4 


2.90 


1.583 


42.4 


1.919 


38.7 


2.152 


-16.8 


3.00 


1.584 


27.6 


1.921 


24.8 


2.150 


-33.5 


3.10 


1.585 


13.3 


1.922 


11.5 


2.147 


-49.8 


3.20 


1.585 


-0.5 


1.922 


-1.4 


2.144 


-65.7 



TABLE V: Model C. 





87r/teo 


= -0.1 


87r/teo = 


87TKeo = 0.1 




M 


uj 2 


M 


J 1 


M 


w 2 


10 15 g cm' 3 


Mq 


10 6 s~ 2 


Mq 


10 6 s- 2 


Mq 


10 6 s- 2 


2.30 


1.471 


77.8 


1.807 


89.5 


2.069 


70.4 


2.40 


1.477 


65.2 


1.816 


77.9 


2.077 


57.0 


2.50 


1.481 


52.8 


1.822 


66.6 


2.082 


43.7 


2.60 


1.484 


40.6 


1.828 


55.5 


2.086 


30.4 


2.70 


1.486 


28.7 


1.832 


44.6 


2.088 


17.2 


2.80 


1.487 


17.1 


1.835 


33.9 


2.089 


4.1 


2.90 


1.488 


5.7 


1.837 


23.4 


2.089 


-9.0 


3.00 


1.488 


-5.5 


1.838 


13.2 


2.088 


-22.0 


3.10 


1.488 


-16.5 


1.838 


3.1 


2.085 


-35.0 


3.20 


1.487 


-27.2 


1.838 


-6.7 


2.083 


-48.0 
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TABLE VI: Model FPS. 







= -0.1 


87TKeo = 


87r/teo = 0.1 




M 


uj 2 


M 


J 1 


M 


w 2 


10 15 g cm" 3 


Mq 


10 6 s- 2 


Mq 


10 6 s~ 2 


Mq 


10 6 s- 2 


2.60 


1.420 


97.3 


1.775 


104.9 


2.039 


68.9 


2.70 


1.425 


81.5 


1.782 


90.5 


2.044 


52.3 


2.80 


1.428 


66.1 


1.787 


76.6 


2.047 


36.0 


2.90 


1.431 


51.2 


1.792 


63.2 


2.049 


20.1 


3.00 


1.433 


36.8 


1.795 


50.1 


2.050 


4.4 


3.10 


1.434 


22.8 


1.797 


37.5 


2.050 


-11.0 


3.20 


1.434 


9.2 


1.799 


25.3 


2.048 


-26.2 


3.30 


1.435 


-3.9 


1.800 


13.4 


2.047 


-41.1 


3.40 


1.434 


-16.6 


1.800 


1.9 


2.044 


-55.9 


3.50 


1.434 


-29.0 


1.800 


-9.3 


2.041 


-70.4 



TABLE IX: Model WS. 





8ttk€o 


= -0.1 


87TKeo = 


87TK£o =0.1 




M 


uj 2 


M 


J 1 


M 


w 2 


10 15 g cm" 3 


Mq 


10 6 s- 2 


Mq 


10 6 s~ 2 


Mq 


10 6 s- 2 


2.30 


1.475 


113.9 


1.812 


118.9 


2.072 


92.3 


2.40 


1.481 


94.6 


1.821 


101.6 


2.080 


74.2 


2.50 


1.486 


76.3 


1.828 


85.3 


2.086 


56.8 


2.60 


1.490 


59.0 


1.834 


69.8 


2.090 


40.2 


2.70 


1.492 


42.5 


1.838 


55.1 


2.093 


24.2 


2.80 


1.494 


26.8 


1.841 


41.1 


2.094 


8.8 


2.90 


1.495 


11.8 


1.843 


27.8 


2.094 


-6.1 


3.00 


1.495 


-2.5 


1.844 


15.0 


2.093 


-20.6 


3.10 


1.494 


-16.2 


1.845 


2.8 


2.091 


-34.6 


3.20 


1.494 


-29.3 


1.845 


-8.9 


2.089 


-48.2 



TABLE VII: Model SLy4. 





S-Kneo 


-0.1 


SuKea — 


87TK60 = 0.1 




M 


uj 2 


M 


w 2 


M 


J 1 


10 15 g cm" 3 


Mq 


10 6 s- 2 


Mq 


10 6 s~ 2 


Mq 


10 6 s- 2 


2.00 


1.692 


143.4 


1.996 


131.5 


2.235 


99.5 


2.10 


1.706 


124.5 


2.011 


113.2 


2.248 


79.1 


2.20 


1.717 


106.6 


2.024 


95.8 


2.257 


59.4 


2.30 


1.725 


89.7 


2.033 


79.2 


2.263 


40.4 


2.40 


1.732 


73.5 


2.040 


63.5 


2.267 


22.1 


2.50 


1.737 


58.2 


2.045 


48.4 


2.268 


4.2 


2.60 


1.740 


43.6 


2.049 


34.1 


2.268 


-13.0 


2.70 


1.743 


29.8 


2.051 


20.4 


2.266 


-29.8 


2.80 


1.744 


16.6 


2.052 


7.4 


2.262 


-46.1 


2.90 


1.744 


4.0 


2.052 


-5.0 


2.258 


-62.0 



TABLE VIII: Model UU. 





8ttk€o 


= -0.1 


8nneo = 


8nne = 0.1 




M 


OJ 2 


M 


w 2 


M 


w 2 


10 15 g cm" 3 


Mq 


10 6 s~ 2 


Mq 


10 6 s~ 2 


Mq 


10 6 s~ 2 


2.00 


1.843 


215.3 


2.127 


173.0 


2.343 


110.2 


2.20 


1.881 


166.5 


2.163 


122.9 


2.366 


51.8 


2.40 


1.905 


124.2 


2.183 


79.3 


2.372 


-0.9 


2.60 


1.920 


87.4 


2.193 


41.2 


2.367 


-48.5 


2.80 


1.928 


54.9 


2.196 


7.8 


2.354 


-91.9 


3.00 


1.931 


26.2 


2.195 


-21.7 


2.338 


-131.5 


3.10 


1.932 


13.0 


2.194 


-35.2 


2.328 


-150.1 


3.20 


1.932 


0.5 


2.191 


-47.9 


2.318 


-168.0 


3.30 


1.932 


-11.3 


2.189 


-60.0 


2.307 


-185.2 


3.40 


1.932 


-22.4 


2.186 


-71.2 


2.296 


-201.8 
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